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Abstract 

We present the results for a search of high-energy muon neutrinos with 
the IceCube detector in coincidence with the Crab nebula flare reported on 
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September 2010 by various experiments. Due to the unusual flaring state of 
the otherwise steady source we performed a prompt analysis of the 79-string 
configuration data to search for neutrinos that might be emitted along with 
the observed 7-rays. We performed two different and complementary data 
selections of neutrino events in the time window of 10 days around the flare. 
One event selection is optimized for discovery of E ~ 2 neutrino spectrum typ- 
ical of 1 st order Fermi acceleration. A similar event selection has also been 
applied to the 40-string data to derive the time-integrated limits to the neu- 
trino emission from the Crab [35]. The other event selection was optimized 
for discovery of neutrino spectra with softer spectral index and TeV energy 
cut-offs as observed for various galactic sources in 7-rays. The 90% CL best 
upper limits on the Crab flux during the 10 day flare are 4.73 x 10 -11 cm -2 
s^ 1 TeV _1 for an E~ 2 neutrino spectrum and 2.50 x 10~ 10 cm -2 s^ 1 TeV _1 for a 
softer neutrino spectra of E~ 2 ' 7 , as indicated by Fermi measurements during 
the flare. IceCube has also set a time-integrated limit on the neutrino emis- 
sion of the Crab using 375.5 days of livetime of the 40-string configuration 
data. This limit is compared to existing models of neutrino production from 
the Crab and its impact on astrophysical parameters is discussed. The most 
optimistic predictions of some models are already rejected by the IceCube 
neutrino telescope with more than 90% CL. 

Keywords: Neutrino and gamma flare, pulsar, nebula. 


1. Introduction 

The Crab supernova remnant, originating from a stellar explosion at a 
distance of 2 kpc recorded in 1054 AD, consists of a central pulsar, a syn- 
chrotron nebula, and a surrounding cloud of expanding thermal ejecta [1]. 
Its bright and steady emission has made it a standard candle for telescope 
calibration. However, the photon emission stability in the X-ray and in the 
7-ray regions is recently being questioned by a number of satellite experi- 
ments. As a matter of fact, a 7% decline of the Crab flux in the 3-100 keV 
region, larger at higher energies, has been observed in the period between 
2008 and 2010 by the Fermi Gamma-ray Burst monitor and confirmed by 
Swift/BAT, RXTE/PCA, and INTEGRAL (IBIS) [2], The pulsed emission 
from RXTE/PCA observations is consistent with the observed pulsar spin- 
down suggesting that the decline is due to changes in the nebula and not in 
the pulsar. 
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The source of energy that powers the Crab is the spin-down luminosity of 
the pulsar. The measured spin-down luminosity of the pulsar is ~ 5 x 10 38 erg 
s' 1 and its rotational period is 33 ms. While a small fraction of this energy 
goes into the pulsed emission, most of it is carried by a highly magnetized 
wind of relativistic plasma, the composition of which is not known. Both pure 

plasma models and a mixture of e± and protons or ions have been proposed 
[1, 3, 4, 6, 10]. The wind terminates in a standing shock and transfers some 
of the energy to accelerating particles. A part of this energy is converted 
into synchrotron emission from radio to MeV y-rays by a population of high 
energy electrons radiating in the nebular magnetic field. The observations 
of the synchrotron emission from the Crab up to the MeV energies, make 
the Crab an undisputed galactic accelerator able to inject electrons up to 
energies ~ 10 15 eV. These high energy electrons inevitably interact with the 
ambient photon fields through inverse Compton scattering, resulting in the 
production of high-energy y-rays observable in the TeV regime [13, 14, 15]. 
The synchrotron emission from the Crab has an integrated luminosity of 
~ 1.3 x 10 38 erg s _1 , that is, at least ~26% of the spin-down luminosity of 
the pulsar is involved in the acceleration of electrons in the energy range 
10 11 - 10 15 eV [1]. On the other hand, the presence of hadrons in the pulsar 
wind and the amount of energy transported by them remain as some of the 
unresolved and interesting questions about the Crab Nebula and plerions in 
general. 

Protons and ions do not lose their energy as efficiently as electrons, and 
hence it is more difficult to observe the products of their interactions. The 
dominant processes, discussed below, are proton-proton and proton-y inter- 
actions, and both processes generate y-rays and neutrinos through meson 
decays. Hence, neutrinos constitute an unique signature for hadron acceler- 
ation while hadronic y-ray production has to be disentangled from inverse 
Compton emission. Hadronic models of the Crab emission assume that the 
pulsar wind is composed of a mixture of electrons and ions. These models 
predict that a significant part of the rotational energy lost by the pulsar 
is transferred through the shock radius to relativistic nuclei in the pulsar 
wind. Relativistic nuclei injected into the nebula can interact with the neb- 
ula matter, and produce cosmic rays and neutrinos via pion decay. Neutrino 
production by protons and nuclei interacting in the pulsar wind in the Crab 
have been discussed in Ref. [3, 4]. According to these models, the nuclei 
can generate Alfven waves just above the pulsar wind shock. These Alfven 
waves will resonantly scatter off and accelerate the positrons and electrons 
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that create the synchrotron emission. In the model described in Ref. [6] neu- 
trinos are produced by heavy nuclei accelerated by the rotating neutron star 
that photo-disintegrate in collisions with soft photons. These models predict 
between 1 — 5 events per year in a cubic-kilometer detector such as IceCube 
when accounting for neutrino oscillations. Inelastic nuclear collisions are con- 
sidered in Ref. [3]. In this paper the predicted rates depend on the Lorentz 
factor, T, of nuclei injected by the pulsar and the effective target density. 
The thermal matter distribution in the Crab is far from being uniform but 
forms filaments. For relativistic protons the effective target density is also 
affected by the structure of the magnetic field in and around these filaments. 
The authors in Ref. [3] provide several expected neutrino fluxes from the 
Crab Nebula as a function of energy, for different assumptions on these two 
parameters. For the highest values of the effective target density, IceCube 
begins to have the sensitivity to probe the highest possible values around 
T < 10' while the favored values of the upstream Lorentz factor of the wind 
are T ~ 10 6 [5], 

Acceleration of positive ions near the surface of a young rotating neutron 
star (< 10 5 yrs) has also been investigated in Ref. [7]. This model describes 
how positive ions can be accelerated to ~ 1 PeV in rapidly-rotating pulsars, 
with typical magnetic fields ( B ~ 10 12 G), by a potential drop across the 
magnetic field lines of the pulsar. Assuming that the star’s magnetic mo- 
ment /i and the angular velocity fl satisfy the relation fl ■ O < 0, protons are 
accelerated away from the stellar surface. Beamed neutrinos (in coincidence 
with the radio beam) are produced by such high energy protons interacting 
with the star’s radiation field when the A production threshold is surpassed. 
Observation of these neutrinos could validate the existence of a hadronic 
component and a strong magnetic field near the stellar surface that acceler- 
ates the charged particles. The predictions in Ref. [8] based on this model 
account for ~ 45 neutrino events/yr from the Crab in a cubic-kilometer de- 
tector in the most optimistic scenario where the fraction of charge depletion 
is assumed to be fd ~ 1/2. In this paper we will show that IceCube data 
severely constrains these optimistic predictions of the model. 

In Ref. [12] a mean prediction of 1.2 neutrino events per year for E v > 
1 TeV was calculated for an underwater cubic-kilometer detector. This pre- 
diction is based on the H.E.S.S. measured y-ray spectrum [13] assuming that 
all the y-rays observed by H.E.S.S. up to 40 TeV are produced by pion decay 
and that the absorption of y-rays is negligible. A similar calculation con- 
necting photon and neutrino fluxes was done in Ref. [9] predicting about 5 
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events from the Crab accounting for neutrino oscillations. For a summary of 
some of the models on neutrino spectra the reader is referred to [10]. 

From Sep. 19 to 22, 2010 the AGILE satellite [16, 17] reported an en- 
hanced 7-ray emission above 100 MeV from the Crab nebula. The flare, 
however, was not detected in X-rays by INTEGRAL [20] observations be- 
tween Sep. 12 and 19 partially overlapping with AGILE observations. It was 
also not confirmed by the SWIFT/BAT [21] in the 15-150 keV range nor by 
RXTE [22] on a dedicated observation of the Crab on Sep. 24. The observa- 
tion was later confirmed by the Large Area Telescope on board of the Fermi 
Gamma-Ray Space Telescope that detected a flare of y-rays (E y > 100 MeV) 
with a duration of ~ 4 days between Sep. 19-22 in the Crab direction [24], 
The observed energy spectrum during the flare interval was consistent with a 
negative power-law with a spectral index of —2.7 ±0.2. The flux increase was 
a factor 5.5 ± 0.8 above the average flux from the Crab. Fermi also detected 
another flare of 16 days in Feb. 2009 corresponding to a flux increase of a 
factor 3.8 ± 0.5 but much softer spectral index (—4.3 ± 0.3). The ARGO- 
YBJ collaboration also issued an ATel on Sep. 2010 on the observation of 
an enhancement of the TeV emission for the same period of time but with a 
wider interval of 10 days. The enhanced TeV emission corresponded to a flux 
about 3-4 times higher than the usual Crab flux in TeV energies [23]. How- 
ever, this observation was not confirmed by MAGIC [25] nor VERITAS [26]; 
Imaging Cherenkov Telescopes in a similar energy range as ARGO-YBJ. The 
spectral and timing properties of the flares indicate that the y-rays are emit- 
ted via synchrotron radiation from PeV electrons from a region smaller than 
1.4 x ICR 2 pc. This dimension is comparable to the jet knots observed close 
to the termination shock of the Crab Nebula [19]. Even though the Crab 
has always been considered to be a source of synchrotron emission, the flare 
represents a challenge to shock diffusive acceleration theory [24], Nonethe- 
less, explanations of the high variability due to electromagnetic phenomena 
have been proposed in Ref. [11] where the emission comes from a part of the 
pulsar wind shock 3 . 

The unusual flaring state of this otherwise steady source, the intensity of 


3 During the final stage of the editing of this paper another large flare was observed 
from the Crab [18]. This flare is even more intense than the one observed in September 
and is being studied by various experiments. Hence, IceCube analysis will happen when 
results from Fermi, other X-ray satellites and other TeV ground based experiments will 
be available. 
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the flare, and the experimental observations in 7-rays motivated this search 
for neutrinos in IceCube in coincidence with the Crab flare of Sep. 2010. 
The IceCube collaboration started a prompt analysis of the then-running 79- 
string configuration. The time window selected for this analysis was the 10 
days interval reported by ARGO-YBJ from September 17 to September 27, 
which contains the Fermi flare window. An unbinned maximum likelihood 
(LLH) method described in Ref. [27] has been applied to search for an ex- 
cess of neutrinos in coincidence with the enhanced 7-ray emission from the 
Crab. The non observation of neutrinos would reinforce pure electromagnetic 
emission scenarios and determine the level at which hadronic phenomena su- 
perimposed on an electromagnetic scenario can be probed. 

The IceCube Neutrino Observatory is a neutrino telescope installed in 
the deep ice at the geographic South Pole. The final configuration comprises 
5,160 photomultipliers (PMTs) [29] along 86 strings instrumented between 
1.5-2. 5 km in the ice. Its design is optimized for the detection of high energy 
astrophysical neutrinos with energies above ~ 100 GeV. The observation 
of cosmic neutrinos will be a direct proof of hadronic particle acceleration 
and will reveal the origins of cosmic rays (CR) and the possible connec- 
tion to shock acceleration in Supernova Remnants (SNR), Active Galactic 
Nuclei (AGN) or Gamma Ray Bursts (GRBs). The IceCube detector uses 
the Antarctic ice as the detection volume where muon neutrino interactions 
produce muons that induce Cherenkov light. The light propagates through 
the transparent medium and can be collected by PMTs housed inside Dig- 
ital Optical Modules (DOMs). The DOMs are spherical, pressure resistant 
glass vessels each containing a 25 cm diameter Hamamatsu photomultiplier 
and its associated electronics. Eight densely instrumented strings equipped 
with higher quantum efficiency DOMs form, together with 12 adjacent Ice- 
Cube strings, the DeepCore array that increases the sensitivity for low energy 
neutrinos down to about 10 GeV. Detector construction finished during the 
austral summer of 2010-11. 

This paper describes in Sec. 2 the data selection, the comparison to simu- 
lation, and the detector effective area and angular resolution for this search; 
in Sec. 3 we summarize the analysis method used; in Sec. 4 the results for the 
flare search are presented. Given the null result, upper limits are provided. 
In Sec. 5 the time-integrated upper limits based on 1 year of data of the 
40-string configuration are presented to summarize what is the impact of the 
IceCube most sensitive limit on existing neutrino production models for the 
Crab. Conclusions are given in Sec. 6. 


2. Data Selection and Comparison to Monte Carlo 

The detection principle of IceCube is based on the charge and time mea- 
surement of the Cherenkov photons induced by relativistic charged particles 
passing through the ice sheet. The PMT signal is digitized with dedicated 
electronics included in the DOMs [31]. A DOM is triggered when the PMT 
voltage crosses a discriminator threshold set at a voltage corresponding to 
about 1/4 photoelectron. Various triggers are used in IceCube. The results 
shown here are based on a simple multiplicity trigger requiring that the sum 
of all triggered DOMs in a rolling time window of 5 /zs is above 8 (SMT8). 
The duration of the trigger is the amount of time that this counter stays at 
or above 8 as the time window keeps moving. Once the trigger condition is 
met, all local coincidence hits are recorded in a readout window of ±10 fis for 
the 40-string run and of 1® /zs (to reduce the noise rate) in the 79-string run. 
IceCube triggers primarily on down-going muons at a rate of about 1.8 kHz 
in the 79-string configuration. Variation in the trigger rate determined by 
atmospheric muons is about ±10% due to seasonal changes [32], Seasonal 
variations in atmospheric neutrino rates are expected to be a maximum of 
±4% for neutrinos originating near the polar regions. Near the equator, at- 
mospheric variations are much smaller and the variation in the number of 
events is expected to be less than ±0.5% [33]. 

For searches of neutrino point sources in the northern sky, IceCube can 
use the Earth as a shield to reduce the background of atmospheric muons 
and detect up-going muons induced by neutrinos. In the northern sky these 
searches are sensitive to neutrinos in the TeV-PeV region. 

In order to reconstruct muon tracks a LLH-based reconstruction is per- 
formed at the South Pole (LI filter) providing a first order background rejec- 
tion of poorly reconstructed events and a selection of high energy muons for 
the southern sky. The data sent through the satellite to the North undergo 
further processing that includes a broader range of more CPU consuming 
reconstructions. This offline processing also provides useful variables for 
background rejection, measurements of the energy and of the angular un- 
certainty, and selects about 35 Hz of the SMT8 data. However, the offline 
processing requires a fair amount of time to be finalized and is not suitable 
for expedited analysis. For the analysis of the Crab flare we used a dedicated 
selection for target of opportunity programs [36]. This online event selection 
and reconstruction is called the online Level 2 filter and selects about 4 Hz 
of data. It provides a reduced data rate (compared to the standard online 
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data) because of stricter cuts than in the offline filter. The loss of sensitivity 
of this stream of data is marginal for E~ 2 neutrino spectra. 

The online L2 filter performs a 8-fold iterative single photoelectron (SPE) 
LLH fit for events with the number of DOMs triggered fewer than 300 and 
a 4-fold iterative SPE fit otherwise. These SPE fits are seeded by a track 
obtained using a single iteration LLH fit [34], While the online Level 2 
selects good quality tracks and high energy muons from the northern sky, 
it is dominated by the background of down-going atmospheric muons and 
therefore further cuts have to be applied before performing neutrino source 
searches. Experimental and simulated data are processed and filtered in the 
same way. The data used for this search concern the period from 2010/08/10 
to 2010/10/12. In this period the detector was running in a stable configura- 
tion. The total live time for that period (considering deadtimes) is 60.9 days. 
Figure 1 shows the data rate of each run included during the selected time 
window as well as the South Pole atmospheric temperature. As can be seen 
at this level, the rate is dominated by down-going atmospheric muons, which 
display larger weather-dependent variations than the final up-going neutrino 
events. 



Figure 1: Data rates per run of the online Level 2 filter of 79 IceCube strings in the time 
window of the Crab flare. The green bar indicates the flaring interval used in this analysis 
according to ARGO measurements [23] . The blue dotted line indicates the temperature 
in the middle stratosphere of the South Pole according to [32]. 

We have performed two dedicated selections starting from the online L2 
filter that we describe below. 
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2.1. Straight Cuts Data Selection 

This dataset is obtained by requiring a good level of reconstruction and 
ensuring degree level accuracy in the tracking errors to reject the misrecon- 
structed down-going atmospheric muons from the real up-going atmospheric 
neutrino sample. The variables used are determined in the offline data pro- 
cessing and have been used for the 40-string point-source analyses in [35] and 
[28]. The final cut level can be achieved by applying the following series of 
cuts on a number of variables to obtain a good agreement between data and 
the simulation of atmospheric neutrinos, with a contamination of the order 
of 5% of atmospheric muons, mainly muons from two cosmic ray showers in 
coincidence in the same readout window. Having these muons with different 
directions gives hit patterns that confuse the reconstruction so that at times 
the result is a misreconstructed up-going track. The cuts are: 

N dlr > 5 ; L dlr > 200 m ; a cr < 5“ ; L„ d < j ™ (1) 

where: 

• N dir : is the number photons detected within -15 and 75 ns with re- 
spect to the expected arrival time of unscattered photons from the 
reconstructed muon-track. Scattering of photons in the ice causes a 
loss of directional information and will delay them with respect to the 
unscattered expectation; 

• L d i r : is the maximum distance in meters between direct photons pro- 
jected along the best muon track solution; 

• cr cr : is the uncertainty on the reconstructed track direction given by 
the LLH-based track reconstruction estimated by a method based on 
the Cramer-Rao inequality [37]; and 

• L re d and L[. ed : are the standard reduced and modified LLH values 
respectively. The reduced LLH is defined as the — log 10 of the LLH 
value of the track reconstruction divided by the number of degrees 
of freedom. The number of degrees of freedom is the number of hit 
DOMs minus five fit parameters, two angles and three coordinates of a 
reference point along the track. It was found by comparing background 
rejection efficiency to signal selection efficiency that a good variable for 
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rejection of background for low energy events is the number of hit 
DOMs minus an effective number of degrees of freedom of 2.5. 

An additional cut to select events in the direction of the Crab (0<>ab = 
122° at the South Pole) has also been applied: ©crab — 10° < 0 rec < 0<>ab + 
10°, where 9 rec is the reconstructed zenith angle of the muon track. No 
further selection in right ascension has been applied. In Tab. 1 the selected 
number of events and the expected number of atmospheric neutrinos and 
muons are given. The final number of events selected for the 10 day window 
of the flare is 354. 

2.2. BDT Data Selection 

The second dataset is obtained by using a multi-variate learning machine. 
In particular this data selection is based on the knowledge and experience 
from previous analyses looking for solar Weakly Interactive Massive Particles 
(WIMPs) with the IceCube detector [38]. During the austral winter the Sun 
is below the horizon at the South Pole and its maximum declination is equal 
to the obliquity of the ecliptic, 23.4°. Since the Crab Nebula lies fairly close to 
the ecliptic plane, the strategies and cuts that are optimized for this specific 
direction can be applied for the Crab direction. 

Starting with the online L2 filtered data selection, as described above, a 
number of additional cuts were applied. The hereby selected events fulfill 
criteria of horizontal tracks passing the detector, to further minimize vertical 
tracks associated with background events. Additionally, the cuts were chosen 
to reduce the tails of distributions of the background into the signal region: 

^travel > “10 777 J (J C OGz < 170 777 ; (J cr < 10° J p av < 150 777 . ; t accu < 3000 HS 

( 2 ) 

where: 

• z travel'- measures the difference in the z positions of the center of gravity 
(COG) of the hits at the beginning of an event (first 1/4 of the hits 
in time) and the COG at the end of the event (last 3/4 of the hits in 
time); 

• <t GO gP is the uncertainty in meters of the ^-coordinate of the COG; 

• p ave : is the mean minimal distance between the LLH track and the hit 
DOMs; and 
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• t accum \ is the accumulation time, defined as the time until 75% of the 
total charge develops in ns. 

Boosted Decision Trees (BDTs) [39], multi- variate learning machines, 
were used in the final analysis step to classify events as signal-like or back- 
ground-like. Eleven event observables, split in two sets of 5 and 6 each, were 
obtained by choosing parameters with low correlation in background (cor- 
relation coefficient |c| < 0.5), but high discriminating power between signal 
and background. The selected observables include N dir , Ldir , cr cr and L' red as 
described within the straight cuts data selection in Sec. 2.1 and z tra vei from 
above. Additionally, observables specifying the geometry, the time evolution 
of the hit pattern, the quality and consistency of the various track recon- 
structions that is defined through the opening angle between the line-fit and 
the LLH tracks, and the number of hit strings are used. Training was done 
with simulated signal events for a soft neutrino spectrum of E~ 3 that also 
well represents the case of an E~ 2 spectrum with a TeV cut-off. A set of 
off-time real data, not used in the flare analysis, was used for training as 
background. The final sample is defined by a cut on the combined output 
(score) of the two BDTs. As in the case of the straight cuts sample, an 
additional requirement of reconstructed zenith tracks within ±10° from the 
Crab has been applied. In Tab. 1 the selected number of events and the 
expected number of atmospheric neutrinos and muons are given. The final 
number of events selected for the 10 day window of the flare is 660 events in 
the northern sky. 


Cut Level 

Data rate 
(Hz) 

Atm. /i rate 
(Hz) 

Atm. rate 

(Hz) 

E- 2 Eff. 
(%) 

Trigger 

1,800 

1,800 

2.59 xlO -2 

- 

Online Level 2 

4.03 

3.11 

7.2 xlO -3 

100 

Straight Cuts 

4.6 x 10” 4 

~ 0 

4.8 xlO -4 

55 

BDT 

8.4 x 10” 4 

~ 0 

8.2 x 10” 4 

61 


Table 1: Data, atmospheric muon, and neutrino expected background rates for different 
cut progression. The signal efficiency for an E ~ 2 neutrino spectrum assuming an emission 
±10° around the Crab with respect to the online Level 2 is also shown. 
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2.3. Comparison Data- Monte Carlo and Detector Performance 

The simulation of atmospheric and signal neutrinos that is used for de- 
termining the selection efficiency, the performance of the detector and to 
calculate upper limits is based on the neutrino generator ANIS [40] and 
the deep inelastic neutrino-nucleon cross sections with CTEQ5 parton dis- 
tribution functions [41]. Neutrino simulation can be weighted for different 
fluxes, accounting for the probability of each event to occur. In this way, 
the same simulation sample can be used to represent atmospheric neutrino 
models such as Bartol [42] and Honda [43] neutrino fluxes from pion and 
kaon decays (conventional flux) and a variety of models for the charm com- 
ponent (prompt flux) [44, 45]. Muons from CR air showers were simulated 
with CORSIKA [46] with the SIBYLL hadronic interaction models [47]. An 
October polar atmosphere, an average case over the year, is used for the 
CORSIKA simulation. Seasonal variations are therefore to be expected less 
than ±10% in event rates [32], Muon propagation through the Earth and 
ice are done using MMC [48]. This simulation is used to verify the level of 
agreement of data and MC from trigger level to Level 1 and to understand 
the level of contamination at final cut level. For the optical properties of 
the ice we used a model obtained from calibrations using the LEDs in the 
DOMs called flashers [49]. This model produces a better agreement between 
data and MC than the model previously used [50]. The simulation propa- 
gates the photon signal to each DOM using light tracking software described 
in [51]. The simulation of the DOMs includes their angular acceptance and 
electronics. The systematic errors on the simulation of the signal used to 
produce the upper limits have been evaluated and presented in Sec. 6 of Ref. 
[35] describing the 40-string time-integrated point source search. The main 
uncertainties on the limits for an E ~ 2 signal of muon neutrinos come from 
photon propagation, absolute DOM efficiency, and uncertainties in the Earth 
density profile and muon energy loss, accounting for a total of 16%. 

Figure 2 shows the data and simulation comparison for some variables at 
the final cut level for the two data samples. As can be seen, the BDT sample 
increases the overall rate by allowing more low quality reconstructed events 
(high L re( f) than the straight cut sample. This is translated into a higher 
neutrino effective area at low energies but also a worse angular resolution as 
can be seen in figure 3. 
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Rates per bin [Hz] Rates per bin [Hz] 


O Data (straight cuts): 4.60e-04 Hz • Data (BDT-wimps): 8.37e-04 Hz 
Atm. v M (straight cuts): 0.00048 Hz Atm. v M (BDT-wimps): 0.00082 Hz 




O Data (straight cuts): 4.60e-04 Hz • Data (BDT-wimps): 8.37e-04 Hz 
Atm. (straight cuts): 0.00048 Hz Atm. (BDT-wimps): 0.00082 Hz 



O Data (straight cuts): 4.60e-04 Hz • Data (BDT-wimps): 8.37e-04 Hz 
Atm. (straight cuts): 0.00048 Hz Atm. v,, (BDT-wimps): 0.00082 Hz 


xIO " 3 



Figure 2: The top-left plot shows the reduced log-likelihood ( L re d ), as defined in section 2, 
distribution for both data (dots) and atmospheric neutrino simulation (green lines) for the 
two data samples. The distribution of the reconstructed energy is shown on the top-right 
plot. The estimated angular error given by the track reconstruction algorithm using the 
Cramer-Rao upper bound is shown on the bottom-left plot while the bottom-right shows 
the azimuth distribution of the final data samples. 
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Neutrino Effective Area (m 2 ) 




Figure 3: Left: Muon neutrino effective area for the two final data samples in a zenith bin 
of ±10° from the direction of the Crab Nebula. Right: Angular resolution defined as the 
median of the point spread function as a function of the neutrino energy for the two data 
samples. The shaded areas represent a ±10% area of the point spread function. 


3. Likelihood analysis 

The method used for this analysis is an unbinned likelihood method [27]. 
This method looks for a localized statistically significant excess of neutrinos 
above the background in the direction of the Crab in coincidence with the 
flare. The same analysis technique has already been applied to AGN flare 
searches in IceCube [28]. The method uses both the reconstructed direction 
of the events as well as an energy proxy, the reconstructed visible muon 
energy, to discriminate any possible signal from background during the time 
interval of the flare. We consider the largest reported time window of 10 days 
by ARGO-YBJ. The applied method describes the data as a two component 
mixture of signal and background. For a data set with N total events the 
probability density of the i th event is given by: 


Us_ 

N 


n. 


S t + (1 - ^)Bi, 


( 3 ) 


where S t is the density distribution for the signal hypothesis and Bi for 
background. The parameter n s is the number of signal events and one of the 
free parameters of the likelihood maximization together with the spectral 
index, 7 , of the signal spectrum distribution. The likelihood of the data is 
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the product of all event probability densities: 


E(n s , 7) 


n 




( 4 ) 


The likelihood is then maximized with respect to n s and 7, giving the best 
fit values n s and 7. The null hypothesis is given by n s = 0 (7 has no meaning 
when no signal is present). The likelihood ratio test-statistic is defined as: 


TS 


-2 log 


- £(n s = 0) - 

- £{n s , %) - ' 


( 5 ) 


The background probability distribution function, or pdf, £>,, is given by: 


Bi = Br c \0 u ^)Br er9y (E h 9 4 ), ( 6 ) 

and is computed using the distribution of data itself. The spatial term 
B^ace^Q .^ .) j g ^] ie even t density per unit solid angle as a function of the 
local coordinates. The energy probability, B^ ner 9 y (Ei, 9 i), is determined from 
the energy proxy distribution of data as a function of the cosine of the zenith 
angle, 6^. This energy proxy, described in detail in [ 35 ], uses the density of 
photons along the muon track due to stochastic energy losses of pair pro- 
duction, bremsstrahlung and photonuclear interactions which dominate over 
ionization losses for muons above 1 TeV. The time probability 6 j) of 

the background can be taken to be flat for this case of a 10 day time interval 
ignoring the seasonal modulations. 

The signal pdf Si is given by: 


5 , = sr ce {\ Xi-x s \,a l )Sr r9y (E l ,6 l , ls )S t r e (U), 


( 7 ) 


where S space depends on the angular uncertainty of the event a* and the 
angular difference between the event position x % from the source position x s . 
The density function S^ ner9V is a function of the reconstructed energy proxy 
Ei , and the spectrum 7^ is calculated from an energy distribution of simulated 
signal in a zenith band that contains the source. The signal time probability, 
Sj tme , depends on the particular signal hypothesis. In this analysis we adopt 
a simple cut in time between t min and t max , which can be expressed as: 


-'time E(t max tij X H(ti 


gtime _ 


t max t min 


( 8 ) 
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where t t is the arrival time of the event, t max and t m i n are the upper and 
lower bounds of the time window defining the flare, and H is the Heavyside 
step function. 

The significance of the result is evaluated by comparing the test-statistic 
with a distribution obtained by performing the same analysis over a set 
of background-only scrambled data sets. The fraction of trials above the 
test-statistic value obtained from data is referred to as the p- value, with 
smaller p- values indicating that the background- only (i.e. null) hypothesis is 
increasingly disfavored compared to the signal-plus-background hypothesis 
as a description of the data. This leads to the definition of the discovery 
potential: the average number of signal events required to achieve a p-value 
less than 2.87xl0 -7 (one-sided 5cr) in 50% of trials. Similarly, the sensitivity 
is defined as the average signal required to obtain, in 90% of trials, a test- 
statistic greater than the median test-statistic of background-only scrambled 
samples. 

4. Results 

The method described in section 3 has been applied to both data samples, 
the one obtained with straight cuts and the one obtained using the BDTs. 
In both cases the best fit resulted in n s = 0 (i.e. an under- fluctuation). 
Figure 4 shows the event distribution for those events with a f > 1, that is, 
only events inside the flare window that contribute to the likelihood. 

As can be seen, due to its higher neutrino efficiency at energies below 
10 TeV the BDT sample has more atmospheric neutrino events. Since the 
background estimation depends on the sample, the signal-to-backgroimd ra- 
tios are different for the same events in the two samples. The highest event 
weight comes from the straight cuts sample. 

Table 2 shows the upper limits set by both data samples for different 
neutrino spectra. Each upper limit is shown both in terms of number of 
signal events that can be rejected at 90% CL, n 9 s 0% , and the flux limit on 
muon neutrinos, for a 9.28 day interval in units of cm -2 s -1 TeV -1 , i.e. 

dN 90% _ <*,90% ( E \~7 
dE V TeV / 

The analysis described and the results given in Tab. 2 rely on the fact that 
background simulation can be performed by scrambling the right ascension 
in real data (even if a signal is present in the data sample the scrambling will 
dilute it over the background). This method of estimating the background 
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Figure 4: Left: Event distribution for the straight cuts for events with ^ > 1. The color 
code represents the event signal-over-background ratio. Right: Same distribution for the 
BDT sample. 

gives robust p - values in terms of systematic uncertainties. Systematic uncer- 
tainties only affect the estimate of the signal flux from the source and the 
upper limits. The systematic uncertainties on the expected flux come from 
photon propagation in ice, absolute DOM sensitivity (±8%), and uncertain- 
ties in the Earth density profile as well as muon energy loss. The main 
uncertainty however is the modeling of Antarctic ice and its effect on the 
photon propagation. In IceCube different ice models have been devised. The 
variation in the upper limits depending on the photon propagation model 
used are within < 10%. Overall the uncertainty on upper limits is 16%. 

5. Impact of IceCube time-integrated limits on models from the 

Crab 

The main goal of the IceCube telescope is the search for cosmic neutrino 
signals that might explain the astrophysical phenomena that give rise to 
the cosmic ray emission. In the absence of detection, constraining models 
can also provide insights about the nature of these phenomena. The best 
available neutrino flux limits for the Crab are based on the time-integrated 
analysis performed during the 375.5 d period corresponding to the 40-string 
configuration of IceCube. We discuss here the impact of these limits on 
different models of neutrino emission from the Crab. Figure 5 summarizes 
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Event weight 


Spectrum 

p — r E cutoff 

(TeV) 

n 90% 

Straight Cuts sample 

3)90% ^min E max 

(GeV) (GeV) 

£90% 

u l* 

n 90% 

<T)90% 

BDT sample 

E min E max 

(GeV) (GeV) 

£90% 

E~ 2 

- 

2.15 

4.84 

IO 33 

1 0 5.7 

1.78 

2.35 

4.80 

10 31 

io 5 - 9 

2.02 

E-2.7 

- 

2.41 

32.6 

10 2 ' 6 

IQ 4 - 9 

6.0 

2.90 

26.3 

IO 2 - 3 

io 4 - 7 

6.77 

E~ 2 

1 

2.80 

309 

io 2 - 4 

io 3 - 5 

21.2 

3.50 

191 

IO 2 - 3 

10 3S 

15.5 

E~ 2 

100 

2.25 

8.59 

10 3 ' 2 

IQS.O 

1.98 

2.51 

8.06 

IO 2 - 9 

io 4 ' 9 

2.05 

E~ 2 

1000 

2.20 

5.52 

io 3 - 3 

io 5 - 6 

1.79 

2.34 

5.31 

IO 30 

io 5 ' 5 

1.86 


Table 2: Upper limits of the Crab Sep. 2010 flare using Neyman for both samples and 
different neutrino spectra including those with an exponential energy cut-off expressed 
as E~ n exp (— E/E cuto ff) where E cu toff is the energy cut-off. The number n 90% is the 
limit in terms of number of signal events for a 90% confidence level and is the flux 
upper limit in units of 10” n cm~ 2 s -1 TeV -1 for a 9.28 days flaring interval. The resulting 
neutrino luminosity limit, is given in units of 10 35 erg s _1 and it was calculated 

by integrating dN 90% /dE x E over the energy range from E min to E max to contain 90% 
signal of the spectrum and multiply by 47 tgP where d is the distance to the Crab Nebula 
(d = 1850 pc). 


a number of different predicted fluxes described in the introduction of this 
paper and where the 40-string configuration limits stand [35]. Upper limits 
are defined as the 90% confidence level (CL) using the method from Feldman 
& Cousins [52], The green line (solid) corresponds to the flux predicted in [12] 
based on the 7 -ray spectrum measured by H.E.S.S. and the corresponding 
upper limit (dashed). The black line represents the estimated flux based on 
the resonant cyclotron absorption model proposed in [3] for the case of a wind 
Lorentz factor of T = 10' and the most optimistic case of the effective target 
density. The red and blue lines represent the two predicted fluxes according 
to [7] for the cases of linear and quadratic proton acceleration respectively. 
The most optimistic version of this model (for both linear and quadratic 
proton acceleration) can be rejected with more than 90% CL using the time 
integrated data from 40 string configuration constraining this way the value 
of the charge depletion fraction. 

6. Conclusions 

Searches for neutrinos in coincidence with the Sep. 2010 Crab flare have 

been presented in this paper. The data used was taken with the 79-string 
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Figure 5: Predicted fluxes and upper limits based on the IceCube 40 string configuration 
on several models from the Crab. Solid lines indicate the predicted flux and dotted lines 
the corresponding upper limit for a 90% CL. The green lines are the predicted flux and 
corresponding upper limit based on the model proposed in [12]. The red and blue lines 
correspond to the model in [7] for the cases of linear (1) and quadratic (2) proton accelera- 
tion. The black line represents the estimated flux for the most optimistic model proposed 
in [3] based on resonant cyclotron absorption model and its corresponding upper limit. 
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configuration of IceCube. This is the first analysis of data taken by this con- 
figuration and represents the first rapid response analysis of IceCube to an 
astronomical event such as the flaring of an otherwise steady standard candle 
source. Two different approaches of event selection have been followed. One 
using direct cuts on quality reconstruction variables and optimized for dis- 
covery for E ~ 2 neutrino spectra, and the other based on multivariate analysis 
and optimized for discovery at lower energies, important for galactic sources 
that have soft spectra with cut-offs at TeV energies. The two data sets how- 
ever showed a background under-fluctuation during the time interval consid- 
ered. The corresponding upper limits based on generic neutrino spectra have 
been shown for the flaring state of the Crab. 

Assuming isotropic emission from the shock (even if this may not be the 
case for a highly relativistic pulsar wind) our limit for E~ 2 corresponds to 
a neutrino luminosity constraint for the flare state of about ~ 2 x 10 35 erg 
s^ 1 , and ~ 1.5 x 10 36 erg s' 1 if a neutrino cut-off of 1 TeV is assumed. 
In both cases the resulting neutrino luminosity constraint is about 2-3 
orders of magnitude lower than the spin-down luminosity of the pulsar and 
comparable to the peak isotropic y-ray luminosity ~ 5 x 10 35 erg s” 1 measured 
by AGILE [17] in the energy range from 0.1 to 10 GeV. 

In addition to the flare analysis we calculated the current best limits set 
by IceCube on different models for neutrino emission from the Crab Nebula. 
These limits are based on the time-integrated analysis of IceCube with the 40- 
string configuration of the detector. The upper regions of the most optimistic 
models can be rejected with more than 90% CL providing useful constraints 
on adjustable parameters of these models. Taking the neutrino spectrum 
derived from the y-ray observations from the Crab, the constraint in neutrino 
luminosity for the steady emission of the Crab is < 1 x 10 35 erg s _1 which 
is a factor ~ 1.7 larger than the luminosity in y-rays assuming the y-ray 
spectrum measured in Ref. [13] integrated over the energy range between 
400 GeV - 40 TeV. 

In the future the IceCube detector will combine datasets from different 
detector configurations. When the different livetimes of the 40-string con- 
figuration data and the full detector will be summed, the sensitivity will 
improve by about a factor of five making this search more predictive. 
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